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Abstract 

An analytical expression for the spin-averaged amplitude squared 
of the rare muon decay e^e~e^V(.v^ is calculated. Monte Carlo 

phase space simulation using the analytical expression for the ampli- 
tude has been used to get various differential distributions of charged 
leptons. The approximate analytical expression for the total energy 
spectrum of charged leptons near the end point is presented. The de- 
pendence of branching ratio on cuts in total energy of charged leptons 
is studied taking into account an experimental energy resolution. It 
is shown that the measured branching ratio is very sensitive to the 
energy resolution. 
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The observation of lepton flavor violation (LFV) processes would indicate 
new physics beyond the Standard Model [1], [2]. The lepton flavor violation 
effects in processes of fi e conversion in a muonic atom, radiative muon 
decay fi^ e'^'-f and muon decay fi^ — > e^e~e^ may be large enough to 
be detected in the future experiments. Ideas [3] of improvement of muon 
beam intensity by a few orders of magnitude lead to new possibilities in 
designing of new lepton flavor violation experiments. Analysis of rare back- 
ground processes is very important for understanding the feasibility of such 
experiments. 

In this article we consider the process 



e+e e+i/eZ/^ (1) 

which is the most important background process in searches for the lepton 
flavor violating muon decay fi^ — > e~^e~e~^. The current measured branching 
ratio for process (1) isT{fi+ e+e-e+iy^i7^)/Ttotai = (3.4±0.4) x 10^^ [4], [5]. 

The detailed analysis of background processes for new experiments search- 
ing for the muon decay /i"*" — > e"'"e~e+ requires knowledge of the amplitude 
of process (1). 

The branching ratio of process (1) was first calculated in [6] and [7] and 
later in [8], [9] and [10]. Unfortunately works [6] and [8] were not published, 
and articles [7] and [9] give only final results of calculations for specific experi- 
mental cuts or without cuts. In [10] the analysis of process (1) was performed 
but the final expression for the squared amplitude was not given. 




Figure 1: Two Feynman diagrams for the process (1). The other two dia- 
grams are obtained by the interchange p ^ p2. 

In order to allow calculation of the branching ratio of process (1) for 
any desired cuts and parameters an analytical expression for spin-averaged 
amplitude squared |Fp is presented in the Appendix. The amplitude was 
calculated for the case of a V-A interaction and unpolarized muons by using 
Mathematica 6.0 and FeynCalc 6.0.0 [11]. The branching ratio calculation 

of process (1) combines |Fp with a phase space Monte Carlo simulation. 
In comparison with the previous articles (such as [9] using the 1.3 x 10^ 
events) the statistics was significantly increased to 10^ and 5 x 10^ events in 
dependence on cuts applied to charged lepton energies. 



In the lowest order in perturbation theory process (1) is described by 
four Feynman diagrams, two of which are shown in Fig.l. The other two 
are obtained by the interchange p ^ p2- Additional diagrams corresponding 
to particle emission from internal W-boson lines are suppressed by a factor 
~ {m^/Mwf ~ 10-6. 

Below we present the results of these calculations of the branching ratio 




Figure 2: The distribution in positron energy dR/dxp (left) and the distribu- 
tion in electron energy dR/dx^ (right) where Xp = Ep/m^ and x^ = Ef,/m^ 




Figure 3: The distribution in positron variable dR/dpdifj (left) and the dis- 
tribution in total energy of charged leptons (right) where y = Etot/fTiii- 

Distributions in positron energy and electron energy without cuts for 
process (1) are shown in Fig. 2. The positron distribution is characterized by 
a long tail up to the maximum positron energy. The electron distribution 
peaks near zero and decreases quickly as electron energy tends to its upper 
limit. 



The distribution in relative positron momentum difference pdiff = \p — 
P2I/I P*^ + P2I where p and p2 are positron momenta is shown in Fig. 3 (left). 
The Pdiff distribution vanishes atpdifj = in agreement with Fermi statistics 
of positrons. A distribution in total energy Etot of all charged leptons without 
cuts is shown in Fig. 3 (right). The spectrum is steep in the region above 
0.5m^. This is important for background suppression in searches for lepton 
flavor violating process fi^ —>■ e~^e~e^. 



Cut (m^ - Etot) 


Statistics 


Branching ratio R 


Irrie 


5 X 10*^ 


(2.83 ±0.16) X 10-19 




5 X 10^ 


(4.660 ±0.046) X lO^^^ 


lOnie 




(3.091 ±0.032) X 10-13 






(7.127 ±0.013) X 10-9 


lOOme 


10^ 


(2.1123 ±0.0022) X 10"^ 


no cut 




(3.5908 ± 0.0033) x 10-^ 



Table 1: The branching ratio R versus cut (m^ — Etot) on the total energy of 
charged leptons. The cut is measured in the electron mass rrie. 



The results of calculations of the branching ratio R are presented in Table 1 
and Fig. 4 (right) versus cut on {m^ — Etot) where Etot is the total energy of all 
charged leptons. Without cuts our result is in agreement with the results of 
articles [7], [9] and [10]: R = (3.54±0.09) x 10-^ R = (3.5916±0.0022) x 10-^ 
and R = (3.6 ±0.1) x 10"^ respectively. For the cut (m^ — Etot) = 5me our 
result is close to those of [9] but for cuts of 10, 50, lOOmg our results are 
less by a factor 10. We note that the differential distribution dR/dEtot of [9] 
contradicts the results quoted in their table. 




Figure 4: The total energy spectrum dR/dy near the endpoint (left) and the 
branching ratio R versus cut on the total energy of charged leptons (right) 
where y = Etot/^n- Smooth lines show the results of fit. 



Also, we would like to note that uncertainties in [9] appear quite small 



for the given number of trials. Our statistics is about 100 times higher but 
uncertainties are comparable with [9]. 

Note that the cut < Smg (see first lines in Table 1) corresponds to the 
region of interest in searches for the lepton flavor violating process /i"*" — > 
e+e~e+. By applying fit to the region near the endpoint (see Fig.4(left)) an 
approximate expression for the spectrum in this region was found to be 

dR/dy = C -{l-yf (3) 

where C = 1.4 x 10"^ and y = Etot/fT^^i- 

Fit of the branching ratio R (see Fig.4(right)) near the endpoint gives an 
approximate expression 

= 2.99 • lO-i'^ (Ihi^Il^Y (4) 
^ me 

Below we show that the steep behavior of the total energy spectrum near 
the endpoint with a finite detector resolution lead to a significant increase 
in the number of background events of type (1) in searches for lepton flavor 
violating muon decay /i"*" — >• e'^e~e'^ . 

The branching ratio to have a signal from muon decay (1) above a thresh- 
old m^ — A is given by 

i?e.p = j dEZ J dR/dEtot ■ f{EZ - Etot)dEtot (5) 



where E^ is the measured total energy of charged particles, Etot is the 
true total energy, A is a value measured from the total energy endpoint, 
dR/dEtot is the total energy spectrum near the endpoint, and f is the reso- 
lution function of the detector. 

For the detector response function of the Gaussian form f{x) — 
^ ^ exp{ — ^ ) the branching ratio Rexp, given by Eq.(5), can be calcu- 
lated analytically by using expression (3) for the total energy spectrum. In 
this case 



Re.p = I- (6) 

6v27r rriij. 



The factor I is given by 



y|(15 + Abu'' + Ibu" + v!')ErJc{^) + ii(33 + Uu" + ^^^)e-"'/^ (7) 



where u = A/a, Erfc{z) = 1 — ^ J dt 



In the limit A » a (an ideal detector) the branching ratio is given by 

■e^P 6 V - 



Rexp — ^{—T- In the case A ~ cr the branching ratio is very sensitive to 



the resolution cr and it is proportional to cr^. Due to the steep total energy 
spectrum near the endpoint the detector resolution changes significantly the 
measured branching ratio Rexp- For example for the typical experimental 2a 
cut on the total energy — Etot\ < 2cr, with cr = A/2 and A = Smg, the 
branching ratio Rexp is increased by a factor 7.8 from 4.66 x 10~^^ (see Table 
1) to 3.64 X 10-^1 
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Appendix. The branching ratio for Muon De- 
cay /i+ e^e~e^Ve'i^^x' 

The kinematical variables of the process (1) were generated by using Phas- 
eSpace.C program of CERN ROOT software package [12] and a weight equal 
to the phase space factor multiplied by the squared amplitude was assign to 
each event. In this program the phase space volume is defined as 

where Pf and are 4-momcnta of final and initial states of the system, 
N is the number of particles in the final state. 

For the process (1) the branching ratio R of the process can be presented 
in the form 



P = Cx|P|2d$ (9) 

where |Pp is the normalized spin-averaged amplitude of the process 
squared. 

The following parameters were taken from Particle Data Group tables [4]: 
muon mass u = 0.105658367 GeV; electron mass m = 0.000510998910 GeV; 
a = 1/137.035999679. 

In terms of these parameters constant C from the previous equation is 
given by 

C = 384Q;V(7rV^). (10) 

Also we define u2 = v?, m2 = m^, mA = m"^, titj = Utj which is the 
scalar product of 4-momenta t(g,p,pi,p2, ^i, ^2) as shown in Fig.l. 
qps = qp*qp; 
qpl2 = qpl*qpl; 
qp22 = qp2*qp2; 
ppl2 = ppl*ppl; 
pp22 = pp2*pp2; 
plp22 = plp2*plp2; 

CI = 1.0/(2.0*(m2 + ppl + pp2 + plp2)); 
C2 = 1.0/(2.0*(m2 - qpl - qp2 + plp2)); 
C3 = 1.0/(2.0*(m2 - qp - qpl + ppl)); 
Dl = 1.0/(2.0*(m2 + plp2)); 
D2 = 1.0/(2.0*(m2 + ppl)); 

trll = -(qk2*(p2kl*(ppl2 - ppl*(m2 + pp2) + m2*(m2 + plp2) 
- pp2*(2.*m2 + plp2)) + plkl*(m4 - m2*pp2 + pp22 + m2*plp2 - 
ppl*(2.*m2 + pp2 + plp2)) + pkl*((2.*m2 - pp2)*(m2 + plp2) - ppl*(m2 
+ 2.*pp2 + plp2)))); 

trl2 = m2*pkl*plk2*qp - m2*plkl*plk2*qp + m2*pkl*p2k2*qp - 
m2*p2kl*p2k2*qp - 2.*m2*pkl*qk2*qp - m2*plkl*qk2*qp 



m2*p2kl*qk2*qp + pkl*plk2*qp*plp2 

+ p2kl*plk2*qp*plp2 + pkl*p2k2*qp*plp2 + plkl*p2k2*qp*plp2 

- 2.*pkl*qk2*qp*plp2 - plkl*qk2*qp*plp2 - p2kl*qk2*qp*plp2 + 
qkl*(m2*qk2*ppl + m2*p2k2*pp2 + m2*qk2*pp2 - p2k2*ppl*plp2 
+ qk2*ppl*plp2 + qk2*pp2*plp2 - 2.*m2*pk2*(m2 + plp2) + 
plk2*(m2*ppl - pp2*plp2)) - m2*pkl*pk2*qpl + m2*plkl*pk2*qpl 
+ pkl*p2k2*ppl*qpl + 2.*p2kl*p2k2*ppl*qpl - p2kl*qk2*ppl*qpl 

- pkl*p2k2*pp2*qpl - 2.*plkl*p2k2*pp2*qpl + 2.*pkl*qk2*pp2*qpl 
+ plkl*qk2*pp2*qpl - pkl*pk2*plp2*qpl - p2kl*pk2*plp2*qpl - 
m2*pkl*pk2*qp2 + m2*p2kl*pk2*qp2 - pkl*plk2*ppl*qp2 

- 2.*p2kl*plk2*ppl*qp2 + 2.*pkl*qk2*ppl*qp2 + p2kl*qk2*ppl*qp2 
+ pkl*plk2*pp2*qp2 + 2.*plkl*plk2*pp2*qp2 - plkl*qk2*pp2*qp2 - 
pkl*pk2*plp2*qp2 - plkl*pk2*plp2*qp2 + klk2*(2.*m2*qp*(m2 + plp2) 
+ pp2*(plp2*qpl - m2*qp2) + ppl*(-(m2*qpl) + plp2*qp2)); 

trl3 = 2.*qk2*(plkl*pp2*(-2.*m2 + pp2) + pkl*(ppl*(m2 - pp2) + 
m2*(m2 + plp2) - pp2*(2.*m2 + plp2)) + p2kl*(ppl*(m2 - pp2) + m2*(m2 
+ plp2) - pp2*(2.*m2 + plp2))); 

trl4 = (m2*pkl*plk2*qp + m2*plkl*plk2*qp + 4.*m2*p2kl*plk2*qp 
m2*pkl*p2k2*qp - m2*plkl*p2k2*qp 

- 2.*m2*pkl*qk2*qp - 2.*m2*plkl*qk2*qp - 4.*m2*p2kl*qk2*qp - 
2.*plkl*plk2*pp2*qp + 2.*plkl*qk2*pp2*qp + 2.*pkl*plk2*qp*plp2 
+ 2.*p2kl*plk2*qp*plp2 - 2.*pkl*qk2*qp*plp2 - 2.*p2kl*qk2*qp*plp2 

- qkl*(-2.*(m2 + ppl)*(m2*p2k2 - qk2*pp2) - plk2*(ppl*(m2 + 
2.*pp2) + m2*(m2 + pp2 - plp2)) + m2*pk2*(m2 + ppl + pp2 + 
plp2)) - m2*pkl*pk2*qpl - m2*plkl*pk2*qpl - 4.*m2*p2kl*pk2*qpl 
+ m2*pkl*p2k2*qpl - m2*plkl*p2k2*qpl + 2.*m2*p2kl*p2k2*qpl + 
2.*m2*pkl*qk2*qpl + 2.*m2*plkl*qk2*qpl + 4.*m2*p2kl*qk2*qpl + 
2.*pkl*p2k2*ppl*qpl + 2.*p2kl*p2k2*ppl*qpl + 2.*plkl*pk2*pp2*qpl - 
2.*p2kl*qk2*pp2*qpl - 2.*pkl*pk2*plp2*qpl - 2.*p2kl*pk2*plp2*qpl + 
m2*pkl*pk2*qp2 + m2*plkl*pk2*qp2 - m2*pkl*plk2*qp2 
+ m2*plkl*plk2*qp2 - 2.*m2*p2kl*plk2*qp2 + 2.*m2*pkl*qk2*qp2 + 
2.*m2*p2kl*qk2*qp2 - 2.*pkl*plk2*ppl*qp2 - 2.*p2kl*plk2*ppl*qp2 + 
2.*pkl*qk2*ppl*qp2 + 2.*p2kl*qk2*ppl*qp2 + klk2*(m2*qp*(m2 + ppl 
+ pp2 + plp2) - (ppl*(m2 + 2.*pp2) + m2*(m2 + pp2 - plp2))*qpl - 
2.*m2*(m2 + ppl)*qp2))/2.0; 

tr22 = -(pkl*(-(plk2*(m2*u2 + plp2*(u2 + qpl) + qpl*(2.*m2 - qp2) 
+ m2*qp2 + qp22)) + qk2*(qpl*(m2 - 2.*qp2) + m2*(m2 + u2 + qp2) + 
plp2*(m2 + u2 + qpl + qp2)) - p2k2*(qpl2 + qpl*(m2 - qp2) + plp2*(u2 
+ qp2) + m2*(u2 + 2.*qp2)))); 

tr23 = (-2.*m2*pkl*plk2*qp + m2*plkl*plk2*qp - m2*p2kl*plk2*qp 
+ m2*plkl*p2k2*qp 
+ m2*p2kl*p2k2*qp + 2.*m2*pkl*qk2*qp + 2.*m2*p2kl*qk2*qp - 
2.*pkl*plk2*qp*plp2 - 2.*p2kl*plk2*qp*plp2 + 2.*pkl*qk2*qp*plp2 + 
2.*p2kl*qk2*qp*plp2 - qkl*(-2.*(m2*pk2 - qk2*pp2)*(m2 + plp2) + 
m2*p2k2*(m2 + ppl + pp2 + plp2) - plk2*(m2*(m2 - ppl + pp2) 



+ (m2 + 2.*pp2)*plp2)) + 2.*m2*pkl*pk2*qpl - m2*plkl*pk2*qpl 
+ m2*p2kl*pk2*qpl - 4.*m2*pkl*p2k2*qpl - m2*plkl*p2k2*qpl - 
m2*p2kl*p2k2*qpl + 4.*m2*pkl*qk2*qpl + 2.*m2*plkl*qk2*qpl + 
2.*m2*p2kl*qk2*qpl - 2.*pkl*p2k2*ppl*qpl - 2.*p2kl*p2k2*ppl*qpl + 
2.*plkl*p2k2*pp2*qpl - 2.*pkl*qk2*pp2*qpl + 2.*pkl*pk2*plp2*qpl 
+ 2.*p2kl*pk2*plp2*qpl - m2*plkl*pk2*qp2 - m2*p2kl*pk2*qp2 + 
4.*m2*pkl*plk2*qp2 + m2*plkl*plk2*qp2 + m2*p2kl*plk2*qp2 - 
4.*m2*pkl*qk2*qp2 - 2.*m2*plkl*qk2*qp2 - 2.*m2*p2kl*qk2*qp2 + 
2.*pkl*plk2*ppl*qp2 + 2.*p2kl*plk2*ppl*qp2 - 2.*pkl*qk2*ppl*qp2 - 
2.*p2kl*qk2*ppl*qp2 - 2.*plkl*plk2*pp2*qp2 + 2.*plkl*qk2*pp2*qp2 
+ klk2*(-2.*m2*qp*(m2 + plp2) - (m2*(m2 - ppl + pp2) + (m2 + 
2.*pp2)*plp2)*qpl + m2*(m2 + ppl + pp2 + plp2)*qp2))/2.0; 
tr24 

= (qpl*(-(m2*p2kl*pk2) - u2*p2kl*pk2 + m2*qkl*pk2 + m2*pkl*plk2 
+ m2*p2kl*plk2 - m2*pkl*p2k2 - u2*pkl*p2k2 + m2*qkl*p2k2 - 
m2*pkl*qk2 - m2*p2kl*qk2 + 2.*p2kl*plk2*ppl - 2.*p2kl*qk2*ppl + 
2.*qkl*plk2*pp2 - 2.*qkl*qk2*pp2 - plkl*(m2*pk2 + m2*p2k2 + 2.*(plk2 

- qk2)*pp2) - 2.*p2kl*plk2*qp + 2.*p2kl*qk2*qp + 2.*pkl*plk2*plp2 - 
2.*pkl*qk2*plp2 + 2.*p2kl*pk2*qpl + 2.*pkl*p2k2*qpl + klk2*(m2*ppl 
+ pp2*(m2 + u2 - 2.*qpl) + m2*(m2 - qp + plp2 - qp2)) - 2.*pkl*plk2*qp2 
+ 2.*pkl*qk2*qp2))/2. + 

u2*((m2*pkl*plk2 

- 2.*m2*pkl*p2k2 + m2*klk2*ppl + 2.*m2*klk2*pp2 - plkl*(m2*pk2 + 
m2*p2k2 + 2.*(2.*plk2 - qk2)*pp2) + m2*klk2*plp2 + 4.*pkl*plk2*plp2 - 
2.*pkl*qk2*plp2 + p2kl*(-2.*qk2*ppl + plk2*(m2 + 4.*ppl) - 2.*pk2*(m2 

- qpl)) + 2.*pkl*p2k2*qpl - 2 *klk2*pp2*qpl)/4.) + 

m2*((2.*m2*qkl*pk2 - u2*qkl*pk2 - 2.*u2*pkl*plk2 + 4.*m2*qkl*plk2 

- 2.*u2*qkl*plk2 - 2.*u2*pkl*p2k2 + 2.*m2*qkl*p2k2 - u2*qkl*p2k2 - 
2.*m2*pkl*qk2 + u2*pkl*qk2 - 2.*m2*plkl*qk2 
4.*m2*qkl*qk2 + 2.*qkl*plk2*ppl + 2.*qkl*p2k2*ppl - 4.*qkl*qk2*ppl + 
2.*plkl*qk2*pp2 - 4.*qkl*qk2*pp2 - 2.*plkl*plk2*qp + 2.*qkl*plk2*qp - 
2.*plkl*p2k2*qp + 2.*qkl*p2k2*qp + 2.*plkl*qk2*qp + 2.*qkl*pk2*plp2 
+ 2.*qkl*plk2*plp2 - 2.*pkl*qk2*plp2 - 4.*qkl*qk2*plp2 + p2kl*(qk2*(- 
2.*m2 + u2 - 2.*ppl + 2.*qp) - 2.*pk2*(u2 - qpl) - 2.*plk2*(u2 - qpl)) + 
2.*pkl*plk2*qpl + 2.*pkl*p2k2*qpl + 4.*qkl*qk2*qpl - 2.*plkl*pk2*qp2 
+ 2.*qkl*pk2*qp2 - 2.*plkl*plk2*qp2 + 2.*qkl*plk2*qp2 
+ 2.*pkl*qk2*qp2 + 2.*plkl*qk2*qp2 + klk2*(-2.*m2*u2 + 2.*pp2*(u2 

- qpl) + 2.*m2*qpl + qp*(2.*m2 + u2 + 2.*plp2 - 2.*qpl - 4.*qp2) + 
2.*m2*qp2 + u2*qp2 + 2.*ppl*qp2 - 2.*qpl*qp2))/4.) + 

u2*m2*((2.*p2kl*pk2 + qkl*pk2 + 3.*pkl*plk2 + 3.*p2kl*plk2 + 
2.*qkl*plk2 + 2.*pkl*p2k2 + qkl*p2k2 - 3.*pkl*qk2 - 3.*p2kl*qk2 - 
plkl*(pk2 + p2k2 + 2.*qk2) + klk2*(6.*m2 + 3.*ppl - qp + 3.*plp2 - 
qp2))/4.); 

tr33 = -(qk2*(plkl*(m4 + m2*ppl - m2*pp2 + pp22 - (2.*m2 + 
ppl + pp2)*plp2) + p2kl*((m2 + ppl)*(2.*m2 - pp2) - (m2 + ppl + 



2.*pp2)*plp2) + pkl*(m2*(m2 + ppl) - (2.*m2 + ppl)*pp2 - (m2 + 

pp2)*plp2 + plp22))); 

tr34 = m2*pkl*p2k2*qp - m2*p2kl*p2k2*qp - plkl*p2k2*ppl*qp 

- p2kl*p2k2*ppl*qp + 2.*plkl*plk2*pp2*qp + p2kl*plk2*pp2*qp - 
plkl*qk2*pp2*qp - 2.*pkl*plk2*qp*plp2 
p2kl*plk2*qp*plp2 + pkl*qk2*qp*plp2 + 2.*p2kl*qk2*qp*plp2 + qkl*(- 
2.*m2*p2k2*(m2 + ppl) + m2*pk2*pp2 + m2*qk2*pp2 + qk2*ppl*pp2 
+ m2*qk2*plp2 - pk2*ppl*plp2 + qk2*ppl*plp2 + plk2*(-(ppl*pp2) + 
m2*plp2)) + m2*plkl*p2k2*qpl - m2*p2kl*p2k2*qpl - pkl*p2k2*ppl*qpl 

- p2kl*p2k2*ppl*qpl - 2.*plkl*pk2*pp2*qpl - p2kl*pk2*pp2*qpl + 
plkl*qk2*pp2*qpl + 2.*p2kl*qk2*pp2*qpl + 2.*pkl*pk2*plp2*qpl + 
p2kl*pk2*plp2*qpl - pkl*qk2*plp2*qpl - m2*pkl*pk2*qp2 
+ m2*p2kl*pk2*qp2 - m2*plkl*plk2*qp2 + m2*p2kl*plk2*qp2 - 
m2*pkl*qk2*qp2 - m2*plkl*qk2*qp2 

- 2.*m2*p2kl*qk2*qp2 + plkl*pk2*ppl*qp2 + p2kl*pk2*ppl*qp2 + 
pkl*plk2*ppl*qp2 + p2kl*plk2*ppl*qp2 - pkl*qk2*ppl*qp2 

- plkl*qk2*ppl*qp2 - 2.*p2kl*qk2*ppl*qp2 + klk2*(plp2*(ppl*qp - 
m2*qpl) + pp2*(-(m2*qp) + ppl*qpl) + 2.*m2*(m2 + ppl)*qp2); 

tr44 = -(p2kl*(-(pk2*(ppl*(u2 + qp) + m2*(u2 + 2.*qp) + (m2 - 
qp)*qpl + qpl2)) - plk2*(m2*u2 + m2*qp + qps + (2.*m2 - qp)*qpl + 
ppl*(u2 + qpl)) + qk2*(m2*(m2 + u2 + qp) + (m2 - 2.*qp)*qpl + ppl*(m2 
+ u2 + qp + qpl)))); 

matr2e = Cl*Cl*Dl*Dl*trll 
Cl*Cl*D2*D2*tr33; 

matr2mu =C2*C2*Dl*Dl*tr22 
C3*C3*D2*D2*tr44; 

matr2emu = Cl*C2*Dl*Dl*trl2 
Cl*C2*Dl*D2*tr23 + Cl*C3*D2*D2*tr34; 

matr2 = matr2e + matr2mu + matr2emu; 

In this notations the normalized spin- averaged amplitude of the process 
squared = matr2. 

An uncertainty of the branching ratio for n trials is calculated by using a 
formula 



Cl*Cl*Dl*D2*trl3 



C2*C3*Dl*D2*tr24 



Cl*C3*Dl*D2*trl4 



+ 



+ 



fW-iWf .1/2 

^= ( — ) 



n 



(11) 



where average weights are given hy W = ^^^-^^Wi/n and W"^ = Z^ILi W^?/^- 



